Statistical Aspects of Ultracold Resonant Scattering 
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Compared to purely atomic collisions, ultracold collisions involving molecules have the potential 
to support a much larger number of Fano-Feshbach resonances due to the huge amount of ro- 
vibrational states available. In order to handle such ultracold atom-molecule collisions, we formulate 
a theory that incorporates the ro-vibrational Fano-Feshbach resonances in a statistical manner while 
treating the physics of the long-range scattering, which is sensitive to such things as hyperfine states, 
collision energy and any applied electromagnetic flelds, exactly within multichannel quantum defect 
theory. Uniting these two techniques, we can assess the influence of highly resonant scattering in the 
threshold regime, and in particular its dependence on the hyperfine state selected for the collision. 
This allows us to explore the onset of Ericson fluctuations in the regime of overlapping resonances, 
which are well-known in nuclear physics but completely unexplored in the ultracold domain. 

PACS numbers: 34.50.-s, 34.50.Cx 
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I. INTRODUCTION 

Resonances have always played a key role in scatter- 
ing experiments across many areas of physics, serving to 
nail down our understanding of the interaction between 
the collision partners. They play an additional role in 
dilute, ultracold atomic and molecular gases, where reso- 
nance positions can be moved relative to the (essentially 
zero) collision energy by means of applied electromag- 
netic fields. This circumstance allows one to control col- 
lision cross sections, as well as mean-field interactions 
in quantum degenerate gases. Dozens of magnetic-field 
Fano-Feshbach resonances have been identified and char- 
acterized in ultracold collisions of various alkali atoms 
[l]; many are now working tools for research in many- 
body quantum physics. In the case of cold collisions of 
alkali atoms, the resonant states differ from the incident 
scattering states by the change of an internal spin. For 
this reason, the number of resonant states remains typi- 
cally small and the resonances themselves usually remain 
well-separated and tractable. 

This situation can be different, however, for collisions 
involving cold molecules, where rotational and vibra- 
tional excitations can also contribute to resonant states. 
Many such resonances have been predicted in theoreti- 
cal treatments of cold molecular scattering pHT2] . While 
the number of resonances naturally grows in this case, 
nevertheless the individual resonances are typically well- 
resolved and manageable in number. This is particu- 
larly evident in cold collisions of molecules with helium 
atoms, relevant to buffer gas cooling, where light masses 
and shallow potential energy surfaces conspire to keep 
the density of resonant states low [I3l [14] . Resonances 
appear to be resolved even in collisions involving light 
objects other than helium such as O2 [2 , Rb+OH [15], 
N+NH [16], or Mg+NH [17 . In relatively "clean" sys- 
tems like these, there remains hope of explicitly identify- 
ing the quantum numbers of resonances, and using them 
to back out accurate potential energy surfaces (PES). 
Indeed, energy resolution afforded at ultra low temper- 



atures may allow for the elucidation of van der Waals 
[T8| IT9] or transition state [20l [21] resonances, impor- 
tant for unraveling chemical reactions when a barrier is 
present. 

There remains, however, a class of heavier molecules 
that have been or will be produced at ultracold temper- 
atures. Notable among these, and the subject of this 
paper, are diatomic species consisting of pairs of alkali 
atoms. When such a molecule collides with another al- 
kali atom, the PES is sufficiently deep that tens of vi- 
brational levels, and hundreds of rotational levels, may 
be energetically accessible. In this case the density of 
resonant states (DOS) may be so high that individual 
resonances may not even be resolved, let alone identified. 
In such a case, it would be worthwhile to understand the 
effect of all these resonances on observed collision cross 
sections. 

Theories relating to high-DOS scattering have long ago 
been formulated, notably in chemistry and in nuclear 
physics. On the one hand, the theory of unimolecular 
dissociation regards the problem in the time domain. If 
a polyatomic molecule is given enough energy to break a 
particular bond, say by absorbing an appropriate photon, 
it does not necessarily immediately dissociate. Rather, 
it can lose energy in many irrelevant degrees of freedom 
until, by accident, sufficient energy lands in the desired 
bond to break it. The theory of this process, known 
as the Rice-Ramsperger-Kassel-Marcus (RRKM) theory 
[22] , expresses the mean rate of dissociation as 
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Here p represents the (very large) density of resonant 
states, while Na represents the (small) number of quan- 
tum states available at the transition state which lead to 
dissociation. 

On the other hand, scattering experiments in nuclear 
theory have inspired statistical ideas of highly resonant 
scattering more in the energy domain. Again, a high den- 
sity of states is expected because of the many strongly 



interacting nucleon degrees of freedom inside a compound 
nucleus. In this case it is typical to treat the energies of 
the resonances (especially if they are individually distin- 
guished) as random numbers with a characteristic mean 
level spacing d = 1/p. Whereas non-interacting energy 
levels are distributed so that their level spacings obey a 
Poisson distribution, instead these strongly interacting 
levels obey a distribution derived by Dyson and Wigner. 
This distribution is regarded as characteristic of spectra 
for systems whose classical analogs are chaotic [23-25 . 

In this energy-domain picture, the resonance widths 
are related, sometimes in a subtle way, to the Hamilto- 
nian matrix elements W^a that couple a bound resonant 
state /i of the collision complex, to a scattering state a 
[26H28] . In the random matrix theory of nuclear scatter- 
ing, these matrix elements are themselves random num- 
bers, distributed about a mean resonance width f . The 
theory identifies two distinct regimes of resonant scatter- 
ing. In the first, T/d <C 1, meaning that the resonances 
are resolvable (though still distributed randomly). In 
the other limit, T/d ^ 1, the resonances overlap. Rather 
than washing out completely, however, the resulting spec- 
trum exhibits "Ericson fiuctuations" on a scale set by 
f itself [^ni EQ]. Both regimes are observed in nuclear 
physics, with Ericson fiuctuations typically appearing at 
higher energies [31]. 

In this article we apply the methods of random ma- 
trix theory to cold collisions within the Wigner threshold 
regime. The object of our study will be atom-diatom col- 
lisions, which possess far fewer degrees of freedom than 
the polyatomic molecules or complex nuclei described 
above. Nevertheless, it has been well-established that the 
same ideas apply to nominally "simpler" systems, even 
to the level of a single electron in a diamagnetic Ryd- 
berg state [3?, "33] or to conductance fiuctuations in a 
semiconductor device [23], in the quantum chaos regime. 

To apply the statistical model to cold collisions, we 
must balance the highly resonant, strongly coupled, lO^K 
energy physics of the complex against the delicate sub- 
mK energy scales of the ultracold. To do this, we exploit 
ideas of multichannel quantum defect theory (MQDT) 
|34^ 35^. This theory makes a clean distinction between 
the physics of the complex, which is pertinent when the 
colliding species are close together; and the physics of the 
long range scattering, which is sensitive to such things as 
the hyperfine states of the atom and molecule, the low 
collision energy, and any applied electromagnetic fields. 
Uniting these two disparate sets of phenomena, we can 
assess the infiuence of highly resonant scattering in the 
threshold regime, and in particular its dependence on the 
hyperfine state selected for the collision. Although the 
multichannel scattering cross sections are derived from a 
fairly realistic framework, we find nevertheless that the 
simple RRKM rate Eq. (IT]) is a useful tool for interpreting 
the results, even at ultracold temperature. 

The present work is outlined as follows. In Section \n\ 
we detail our theoretical framework, which is divided into 



ing framework and the treatment of the long-range inter- 
actions via MQDT. In Section |IIB we then present our 
approach of treating the highly resonant short-range part 
by means of a statistical framework derived from random 
matrix theory. The essential input parameter for the sta- 
tistical theory is the density of states for the short-range 
resonances; in Section |IIC| we provide estimates for all 
non-reactive A + AB alkali dimer pairs. The question of 
including the density of states due to the nuclear spin de- 



grees of freedom is addressed in Sec. |IID| In Section [III A 
we present exemplary elastic cross sections within the 
Wigner threshold law regime that are derived from our 
theoretical framework. Two particular examples are cho- 
sen: K + LiK, where resonances remain well-separated, 
and Rb + KRb, where the DOS is high. Also, mag- 
netic field dependent thermal rates are provided. Section 
III B| finally shows that ultracold atom-molecule collisions 
demonstrate the onset of Ericson fiuctuations on a com- 
pletely different energy scale than in nuclear physics. In 



two aspects. Section II A introduces the general scatter- 



Section IV we comment on what might be learned from 
experimental data by comparing to the predictions and 
assumptions of our model. With Section [Vj we provided 
a brief conclusion and an outlook on further directions 
for our theory of highly resonant scattering. 



II. THEORETICAL FRAMEWORK 

We deal here with the three-body physics of ultracold 
alkali atoms, a calculation that could, in principle, be 
performed in substantial detail [36l[37]. It is, however, 
an immense labor, and the results, while qualitatively 
meaningful, are unlikely to be quantitatively accurate. 
Even in cases where the calculations can be converged, 
the relevant potential energy surfaces are not known to 
sufficiently high accuracy for ultracold collisions. Nev- 
ertheless, the scattering framework is standard. In this 
section we develop this framework, including our approx- 
imate, statistical version of the resonant states. 



A. Scattering frameAvork 

We begin with a diatomic molecule AB (where A and B 
denote alkali atoms) in its ^S electronic ground state, its 
V = vibrational ground state, and its n = rotational 
ground state, according to the Hund's case b) coupling 
scheme. Examples of such molecules have been produced 
in gases of order /iK temperatures [38H4T] . The molecules 
may or may not be also prepared in their ground state 
of nuclear spin / [42] . These molecules will collide with 
another alkali atom C (typically one of A or B), in its 
'^S ground state, and with its own accessible hyperfine 
degrees of freedom. 

For the time being, we consider cases where chemical 
reactions are not energetically allowed at ultralow tem- 
perature. Therefore, the only collisions we consider are 
those that can change the nuclear spin quantum num- 



bers. Generally, we are interested in the regime where 
the atom's spin state can be labeled by |/, m/), even in 
the presence of a magnetic field. For the molecule, we as- 
sume a magnetic field sufficiently large that the nuclear 
spins Ia and Ib are decoupled and the states can be char- 
acterized by their individual projections on the magnetic 
field axis, |/^My^, /^M^). The observables then consist 
of the collision rate constants 



K 
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where v is the relative velocity before the collision, and a 
is the collision cross section. We omit the nuclear spins 
I A and Ib in Eq. Q since they are not subject to change 
in the collisions we are considering. Assuming a decom- 
position into partial waves \LMl) for the relative motion, 
the cross section is given by 
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in terms of the scattering matrix elements Sa'a- The S- 
matrix describes the possible re- arrangement of angular 
momentum during the collision, but must conserve the 
total projection, M = M'^^ ^ M'^ ^ m'j ^ M'^ ^ m'^ = 
Ma + Mb + mf + Ml + rrin^ with the quantization axis 
applied along the magnetic field direction, if any. In the 
above equality we included the projection m^ of the ro- 
tational quantum number n of the molecule, which is 
needed when considering possible resonant states. For 
the incident and outgoing channels, however, we will al- 
ways assume the ro- vibrational ground state, i.e., v = 
n = rUn = 0. For notational convenience we hereafter 
denote these scattering channel indices as 
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Calculation of a schematic but realistic Sa'a^ including 
its energy- and magnetic field-dependent resonance struc- 
ture, is the goal of this article. 

To achieve this goal, we exploit the conceptual differ- 
ence between the spin channels \a) that describe physics 
at large interparticle separation R] and the numerous res- 
onant states l/i) that differ by rotational and vibrational 
quantum numbers from a, and that describe states of 
the scattering complex. This general separation of states 
is illustrated schematically in Fig. [l] For separations R 
greater than some characteristic distance Rm^ the chan- 
nels a are assumed to be independent of one another, 
and described by simplified long-range interactions of the 
form 
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where Ea{B) is the threshold of the ath channel, which 
may depend on a magnetic field B. Here, rrir is the re- 
duced mass of the scattering partners and Cq is their van 




FIG. 1. (color online) Schematics of our MQDT approach (not 
to scale). In the long-range part, R > Rm, we only consider 
the ro- vibrational ground state of the molecule, but include all 
Na atomic and molecular hyperfine states and different partial 
waves (not shown). In the asymptotic region, R ^ oo, No of 
them are energetically open and Nc are closed. The MQDT 
treatment transforms the short-range i^-matrix K^^ , defined 
at Rm, into a physical scattering matrix S_^^y^ from which 
quantities such as elastic and inelastic cross sections can be 
deduced. K^^ includes the information on the ro- vibrational 
resonances which occur in in the short-range part, R < Rm- 



der Waals coefficient, which is taken to be isotropic in 
this model. 

Dividing the scattering process into short- and long- 
range parts forms the basis of quantum defect theory 
(QDT). Here, we utilize a multichannel formulation of 
QDT along the lines of Ref. [35 . The key feature of 
MQDT is that - once the MQDT parameters have been 
determined for a given class of potentials - one only needs 
to provide the reactance matrix K^^ = i{l — S_^^){1 + 
5^^^)~^ which is deffned at the matching radius Rm be- 
tween the short- and the long-range. The MQDT for- 
malism as outlined in this section then takes care of the 
propagation for R > R^ and directly yields the physical 
scattering matrix S_^ ^^^ which defines the solution vec- 
tors i/?^'^^ of the coupled channel equations for the whole 
scattering process. 
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f^ = ^y2mr/7Th?kae^'^^^''■^ ^a7r/2) ^^^ outgoing (+) and 
incoming (— ) spherical waves, respectively. Once deter- 
mined, S_^^^^ can easily be converted into various observ- 
ables describing the scattering process. 

We apply this formalism explicitly only to the small 
number Na of hyperfine channels belonging to the ro- 
vibrational ground state of the molecule and the ground 
electronic state of the atom. Of these, some number No 
will be energetically open, meaning that for these chan- 
nels E > Ea and the collision partners can escape to 
infinity. The remaining Nc = Na — No closed channels 
do not contribute directly to the physical scattering ma- 
trix, and must be "eliminated" by the usual algebraic 
procedures of MQDT. 



To do so, the short-range if-matrix K^ 
into its open and closed channels at Rm^ 






is partitioned 



(7) 



The closed channels are eliminated in the MQDT sense 
through 
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(8) 



where /3 is a closed-channel MQDT parameter [35 . The 
modified reactance matrix K^ has dimension No x No and 
shows the potential influence of closed channel pathways. 
The transformation to an energy-normalized, nonana- 
lytic long-range representation is achieved by 



K = A^K{l^QK)~^A^- 



(9) 



The physical scattering matrix is finally formed by 

A, Q_^ T]^ and P are diagonal matrices in the asymptotic 
channel space, consisting of the relevant MQDT param- 
eters. The latter are determined as in Refs. [35l[43]. In 
the present from of MQDT we encounter two sets of long- 
range reference functions: (/^,^^) are smooth, analytic 
functions of energy whereas (/, g) are energy-normalized 
but non-analytic functions of energy. (/, g) are solutions 
of the Schrodinger equation in presence of a long-range 
potential V^^{R) and are related to the energy- normalized 
spherical Bessel and Neumann functions via 
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(11) 



Equation (11) defines the MQDT parameter 77. The pa- 



rameter /3 is a negative energy phase that represents the 
phase accumulated in V^^{R). The energy- normalized 
base pair (/, ^) is related to energy-analytic base pair 
{f^->9^) through the transformation 
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(12) 



which defines the MQDT parameters A and Q. 

MQDT has been a hugely successful tool for organiz- 
ing apparently complex spectra of atoms [44 and in sim- 
ply describing resonant scattering, both at thermal ener- 
gies dH [46] and in ultracold atom collisions [35l I47H49] . 
Much of its appeal in these circumstances lies in the fact 
that a matching radius Rm can be chosen so the channels 
that will be closed as i? ^ 00 remain classically open at 
R = Rjn. If this is so, the short-range K matrix be- 
comes a weakly energy-dependent quantity, and complex 
spectra can be unified by the simple algebraic procedures 
described above. For molecular scattering, it remains to 



be seen whether this same simplicity occurs, since for any 
Rm there may be many channels that are already clas- 
sically closed, and hence, impart resonant structure to 
K^"^ . Indeed, there is already some hint in applications 
of MQDT to cold molecule collisions that Rm must be 
chosen carefully to maximize the simplicity of K^^ [12] . 

In the present case of highly resonant scattering, we 
in fact approach quite the opposite limit, where for any 
reasonable Rm most of the resonant ro- vibrational chan- 
nels are already closed. Thus, our IC''^ will necessarily 
be highly energy dependent, exhibiting already the reso- 
nances of interest. Although it is difficult to compute, it 
remains nevertheless a well-defined quantity in the the- 
ory. For our present purposes, we employ MQDT as 
quick, algebraic solution to producing scattering matrices 
S^ ys fQj. Q^ given IC'^ . Arriving at a physically reasonable 
K^"^ is the task we turn to next. 



B. Statistical short-range K-Matrix 

In treating the long-range collision physics by means 
of MQDT, the only quantity left to be determined is the 
short-range i^-matrix. It is indexed by the Na asymp- 
totic channels a, but is influenced by the myriad (i.e., 
N ^ No) of resonant states ji. Quite generally, it can be 
expressed as [31] 



N 
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E - E, 



(13) 



Equation (13) is expressed in the eigenspace of a short- 
range Hamiltonian IP"^ that gives rise to the (unper- 
turbed) short-range levels at E^. Wa^ = W^a ^ire 
(assumed-to-be energy independent) coupling matrix el- 
ements between resonance fi and asymptotic channel a. 
The mean coupling strength of the ath asymptotic chan- 
nel to the short-range resonances is given by the dimen- 
sionless parameter 



R'S^ = IpTa, 



(14) 



^N 

^/i= 



liy^aP is the zero-order aver- 



where f , = (2^/7V) J^^^^ 
age partial width to the decay channel a [28 and p is the 
DOS of the resonances, evaluated at the incident energy. 
The input parameters for the resonant scattering the- 
ory, Eq. (13), are the zero-order positions E^ of the reso- 
nances and the coupling elements Wa^ to the asymptotic 
channels. However, both are usually unknown unless the 
short-range part is known with high precision. To provide 
Eq. (13) with reasonable input parameters, we utilize a 
statistical model where E^ and Wa^ are taken as random 
variables. This model follows closely the random matrix 
theory approach in nuclear reaction physics [31 and can 
be also found in theoretical works on quantum transport 
[23] . as well as in the theory of chemical reactions [28]. 
By employing such a model, we assume that the collision 



complex corresponds classically to a long, chaotic tra- 
jectory that ergodically explores a large portion of the 
allowed phase space. 

Acknowledging the statistical nature of the short-range 
resonance levels, we apply random matrix theory to K^^ 
based on the Gaussian Orthogonal Ensemble (GOE) ac- 
cording to [31 . In particular, we assume that the spec- 
trum E^ of the resonant states is determined by a Hamil- 
tonian H^^^ that is a member of the GOE. As a such, 
the nearest neighbor distribution of the spectrum satis- 
fies the Dyson- Wigner distribution 



^(V) = P(s) 
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(15) 



where 5^ = |£^^+i — E^\/d is the nearest-neighbor level 
spacing in units of the mean level spacing [24 . In 
practice, we produce the spectrum E^ for a given DOS 
p = 1/dhy constructing first a set {s^} of nearest neigh- 
bor splittings satisfying Eq. ([l5| ^; the spectrum is 
then given by E^ = Eq ^ XlitTi "^m? where Eq is an ap- 
propriately chosen offset. An exemplary GOE spectrum 
for Rb + KRb scattering is reproduced in Fig. [2] along 
with its nearest neighbor distribution. 

Since K^^ is expressed in the frame where H^^^ is 
diagonal, the coupling matrix W_ becomes a random pro- 
cess itself [31]. More precisely, its elements are given 
by uncorrelated, Gaussian-distributed random variables 
with mean and variance v'^. Hence, f ^ = 27ri^^ are the 
mean zero-order partial widths. From Eq. (14) we find 
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(16) 



Thus, in order to describe the short-range physics within 
the statistical model, it is sufficient to specify the DOS 
p of the short-range resonances and the mean coupling 
strength Ra to the asymptotic channels. In the present 
work, we will usually assume Ra = 1 for which the 
transmission coefficient T^ between the short and long- 
range channels [28^ 22l EI] , 



Ta 
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(17) 



reaches unity. In other words, outbound flux that has 
left the collision complex and reaches i?^, is assumed 
to continue out with unit probability. Some of this flux 
will later be reflected back to small R due to details of 
the hyperfine channels \a). This effect, however, is fully 
accounted for in MQDT. 

Having a unit transmission probability corresponds to 
the RRKM limit of transition state theory, reached for 
barrier less reactions. In transition state theory, the de- 
cay rate of a metastable state (here: the short-range res- 
onances) is proportional to the ratio between the number 
Na of open channels (here: the asymptotic channels) and 
the level density p of the metastable states ^22^. Indeed, 




v=2,n=0,l, 



v=l,n=0,l,... 
v=0,n=l,2,... 



incoming / outgoing channel: 

rovibrational ground state 

v=n=0 



FIG. 2. (color online) Schematic overview (not to scale) of 
the origin of the short-range resonances and their distribu- 
tion, (bottom) The atom-molecule potential is modeled by 
a Lennard- Jones potential, Eq. (19). The resonant channels 
stem from ro-vibrationally excited states of the molecule. For 



ultracold temperatures, the incident and outgoing scattering 
channels are restricted to the ground ro-vibrational state of 
the molecule. Within this ground state, the various spin states 
are treated explicitly by means of MQDT. (top) Exemplary 
short-range spectrum for s-wave collisions of Rb + KRb; it is 
constructed to satisfy the Dys on- Wigner distribution for the 
nearest neighbors, cf. Eq. (15). 



one recovers for the decay rate /c. 



k = 






2iThp 2'Khp 



= h 



RRKM 



?(0) 



(18) 



when Ta = 1 [28]. In general, however, Ra can act as a 
fitting parameter to real spectra, revealing further details 
on the short-range physics. 



C. Ro-vibrational density of states 

Having the statistical model for the short-range K- 
matrix in hand, Eq. (13), the question of computing the 



mean density of states itself remains. Since in general we 
do not know the short-range potential in detail, it is im- 
possible to calculate the real atom-diatom ro-vibrational 
spectrum, which would give rise to the short-range reso- 
nances in question. We thus pursue the following strat- 
egy to get an adequate estimate of the short-range ro- 
vibrational density of states prv (see also Fig. [2]): 

1. The short-range interaction is approximated by a 
Lennard- Jones potential along the reaction coordi- 



nate R and add the centrifugal barrier due to the 
end-over-end rotation angular momentum L of the 
atom and the molecule about one another, 



V^^\R) 
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(19) 



A fairly realistic estimation of the long-range be- 
havior of this potential is to assume C6,b+ab = 
C6,AB + C'6,B2 [51 J a convenient compilation of the 
Cq coefficients of all alkali dimers can be found in 
[5? . and we use these here. The C12 coefficient in 



Eq. (19) refers to the short-range behavior of the 



potential. It can be expressed in terms of the over- 
all depth De of the potential and the Cq coefficient 
via C12 = C|/4I)e; in the present work, we em- 
ploy the realistic depths De as calculated in [53 . 
De refers to the dissociation energy of the ground 
state trimer AB2 into AB + B. 

2. For every partial wave L of interest, we calculate 
the bound state energies Ea of Fsr (^)- ^ labels 
the vibrational quantum number in i?, and each of 
these states represents a possible short-range reso- 
nance. 

3. Each asymptotic ro- vibrational channel (v, n) can 
give rise to such short-range resonance states. With 
each of these channels we associate a set of reso- 
nance energies Ea '^'"^ that is offset by the corre- 
sponding channel threshold, i.e., Ea '^'^ = Ea + 
Ey^n- The total angular momentum J = L + n is 
assumed to be conserved in the usual quantum me- 
chanical way, hence the triangular conditions hold 
for the possible combinations of L and n. In par- 
ticular, for J = (s-wave collisions), only L = n is 
possible. 

4. All allowed energies Ea '^ form a total spectrum. 
From this spectrum we can extract the mean level 
spacing d and the level density prv = ^/d associated 
with the ro-vibrational resonant states. In doing 
so, we restrict ourselves to a certain energy interval 
centered around the incoming channel threshold. 

5. The analysis so far does not account for any de- 
generacy of the energy levels. We consider the case 
where only the total magnetic quantum number M 
is conserved. Since M = Mx^M^^mf^ML^rrin^ 
there are numerous possibilities to couple to a given 
total M. These degeneracies are accounted for 
when calculating the final DOS. 

Specific examples of DOS calculated in this way are 
presented in Tables [l| and [ill The former shows the depen- 
dence of the resulting DOS on the maximal ro-vibrational 
quantum numbers '^max^^max used in the estimate; the 
Rb -h KRb collision is chosen as a particular exam- 
ple. The DOS increases with the number of allowed ro- 



vibrational levels V. 



and saturates if a large num- 



max5 ^max 

ber of ro-vibrational levels is included. This saturation is 



TABLE L Ro-vibrational DOS prv(mK-^) for ^^Rb+^°K^^Rb 
collisions as a function of the maximal allowed vibrational 
(rows) and rotational (columns) quantum numbers. The to- 
tal angular momentum J = L + n is assumed to be conserved; 
the values reported here are for J = 0, Ml -\-mn — 0. Includ- 
ing more than 25 vibrational and 100 rotational level does not 
increase the DOS further. 
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0.17 
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2.53 
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0.18 


0.46 


1.20 


2.22 


3.46 


4.81 
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0.49 


1.29 


2.34 
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TABLE II. Properties characterizing the atom-molecule po- 
tential, Eq. (19), and the resulting ro-vibrational DOS for 
various atom-molecule pairs. The DOS are calculated for 
L = 0. We picked the isotopes ^Li, ^^Na, ^°K, ^'^Rb, and 
^^^Cs. The X^i;+ molecular PES needed for the calculation 
of the ro-vibrational states of the AB molecule are taken from 
the references provided in the last column. 
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Na+LiNa 


2872 


3015 


0.19 


0.569 


0.11 


m 


K+LiK 


2124 


6190 


0.33 


0.184 


0.06 


[55] 


Rb+LiRb 


2036 


7159 


0.76 


0.056 


0.04 


|56| 


Cs+LiCs 


2502 


9785 


1.41 


0.026 


0.04 
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K+NaK 


1851 


6297 


1.22 


0.150 


0.18 


m 


Rb+NaRb 


1752 


7273 


2.73 


0.050 


0.14 


|59| 


Cs+NaCs 


2186 


9910 


4.96 


0.024 


0.12 


m 


Rb+KRb 


1684 


8798 


4.95 


0.041 


0.20 


ISB 


Cs+KCs 


1821 


11752 


7.68 


0.020 


0.16 


|62| 


Cs+RbCs 


1825 


12135 


12.11 


0.017 


0.21 


m 



because of the finite sampling interval for the calculation 
of the mean DOS: the bound states belonging to highly 
excited ro-vibrational molecular levels simply lie outside 
the sampling region. We remark that the final densities 
of states are rather insensitive to the particular sampling 
interval chosen; in the present work, we use an interval 
of ±5 K centered at the ro-vibrational ground state. In 
addition, as can be seen in Table [T| the dependence of 
Prv on 'i^max^^max is rather weak. Thus, for example, if 
atoms in the collision complex can for some reason only 
access states up to nmax = 60 rather than nmax = 100, 
this would only change our estimate by a factor of two. 
In general, we expect something like one ro-vibrational 
resonance per mK for 5-wave scattering of Rb + KRb. 
This estimate assumes that the molecular states of AB 
are all in their singlet electronic manifold. We estimate 
that including the (much shallower) triplet states would 
increase the DOS by ~ 10%. Therefore, we do not con- 
sider these states in the present work. 



We provide similar estimates of ro- vibrational DOS for 
various collision partners in Table [llj assuming that all 
energetically allowed v and n states contribute. We also 
include in this table the basic molecular data from which 
the DOS estimates were obtained. As one might expect, 
Prv is larger for heavier collision systems. In particu- 
lar, the DOS for Cs + RbCs collisions is two orders of 
magnitude higher than for the light Na + LiNa collision 
complex. A useful way to express prv is in units of states 
per van der Waals energy £^vdw = ^^(2m^)~^/^Cg , 
see also Table [lIj In this representation, the larger van 
der Waals energy scale of lighter molecules compen- 
sates for their smaller number of bound states. In the 
end, all considered DOS are roughly the same, namely, 
Prv(£^7dw) ^0-1 within a factor of 3. 

The DOS provided in Tables [l| and |n| are specific ex- 
amples for J = 0, for which L = n needs to be satis- 
fied in order to conserve the total angular momentum 
J = L + n. For J y=^ 0, there are 2 J+ 1 possibilities for L 
and n to couple (L = n, L = n±l, ..., L = n ^ J) and 
therefore the DOS increases by approximately the same 
factor. 

We remark that in the above considerations the pres- 
ence of nuclear spin degrees of freedom is not taken into 
account. In the following subsection we will therefore 
discuss the possible infiuence of the spin on the densities 
of states. 



D. The role of nuclear spins 

Thus far we have considered only the density of states 
due to rotations and vibrations, i.e., due to the relative 
motion of the three alkali atoms, denoted by p^^. The 
DOS will multiply, however, if the nuclear spin degrees 
of freedom become involved. To see whether the nuclear 
spin may change during the collision, we employ a semi- 
classical analysis as follows. Once the collision complex 
is formed, it lives, on average, for an amount of time r 
that is related to the mean resonance width by r = h/T. 
During this time, the nucleus of any given atom follows 
a chaotic trajectory through phase space, according to 
our ergodic assumption. The nuclear spins are infiuenced 
during this time by a hyperfine Hamiltonian JT^hf that 
varies in time as the collision complex explores the phase 
space. The dominant part of this Hamiltonian arises from 
the magnetic dipole interaction of the nuclear spin with 
the spin of the electron immediately in orbit above it in 
the same atom. The electron spin, however, is subject 
to fluctuations during this classical trajectory. We there- 
fore expect that the nuclear spin experiences a rapidly 
time- varying change in its Hamiltonian, SHiif{t). 

The magnitude of these fluctuations can be estimated 
by the following argument. We regard the collision com- 
plex semiclassically as a repeated set of mini-collisions, 
occurring at average time intervals At. For instance, at 
one moment the complex might resemble A+ (BC)*, i.e., 
the A atom is loosely bound to a molecule BC, which 



is excited into some ro-vibrational state. Because BC is 
excited, A cannot escape, but rather returns to collide 
again. This collision might result in a different complex, 
say B + (AC)*. Consider then the nucleus attached to 
atom A. Before this collision, this nucleus experiences the 
unpaired electron on the A atom, and hence essentially 
the entire hyperflne interaction determined by the cor- 
responding magnetic dipole constant Ahf [64 . After the 
collision, the atom A is locked into a singlet state with 
atom C, and the nucleus sees no hyperflne interaction 
at all, apart from the modest nuclear quadrupole inter- 
action. Thus, the hyperflne interaction experienced by 
any given nucleus in the complex is effectively switched 
randomly between full strength and zero, at random in- 
tervals ~ At. 

Let Jiwi2 denote the energy difference between two nu- 
clear spin states in the absence of these fluctuations. 
Then a nucleus initially in one of the states will end up 
in the other at time r with a probability amplitude 
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in terms of the Fourier transform ^i^hf(^i2)- The per- 
turbing Hamiltonian SHiif{t) will fluctuate on a charac- 
teristic time scale At, set roughly by the mean collision 
time of an atom in the complex with another atom. The 
Fourier transform of ^i7hf(t) is then nonzero only over 
some flnite bandwidth Q = 2i\ j At. To conserve the in- 
tensity of the fluctuations over r in the time domain and 
Vt in the frequency domain, the root-mean-squared av- 
erages of the fluctuation and its Fourier transform must 
satisfy 



v/([.5i/hf(t)]2)tv^ « ^([5iIhfH]2)^v^. (21) 

If we assume that the time domain fluctuations are ran- 
dom white noise, then the power spectrum is approxi- 
mately independent of frequency within the bandwidth 

1], i.e., ^^hf (c^) = ^iiti = const, and J {[6Hm{uj)]^) = 
^^hf correspondingly. In particular, we assign it 
this value at the transition frequency, (^^hf(<^i2) = 
SH^^K Employing c{t) = V27r /ih x ^^hf(^i2) and 
y^{[SHii^{t)]'^)t = Ahf/2, we flnd for an estimate of the 
transition probability 
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The flrst factor in parentheses denotes the size of the per- 
turbation, times the length of time is acts, which would 
be the probability that a smoothly-varying perturbation 
changes the spin state. The second factor in parentheses 
accounts for the fluctuations on a time scale At. Once the 
perturbing Hamiltonian takes a certain value, the nuclear 



TABLE III. Mean collision time At [Eq. J24|)], lifetime r, 
nuclear spin transition probability P [Eq. (|22[)j, and nuclear 
spin enhancement factor A^nuc [Eq. (26)] for various collision 
complexes in case of L = 0. 
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spin has only a time At to respond to this perturbation, 
i.e., by precessing around the instantaneous local mag- 
netic field. After time At, the perturbation randomly 
switches to something else, and the nuclear spin attempts 
to follow a new local field. If At is much smaller than 
the nuclear-spin-changing period h/uji2 (as it is in our 
case), then the nuclear spin has a hard time changing at 
all; more rapid collisions actually reduce the transition 
probability. On the other hand, if the collisions occur 
rarely on the time scale h/uji2^ then the relevant At is 
reciprocal to the hyperfine interaction itself. In this case, 
the second factor in (22) is unity, and we reduce to the 



familiar case of a slowly-varying perturbation. 

Representative values of P in case of s-wave collisions 



are given in Table III for some of the atom- molecule pairs 
we are considering. For calculating these values, we esti- 
mate the complex lifetime by means of its RRKM value, 
r = h/T = 27rfiprv, where prv is the ro- vibrational DOS 
of the collision complex as provided in Table [ll| for L = 0. 
The mean collision time can be estimated by averaging 
over classical trajectories in a pure Cq potential as fol- 
lows. The time the atom needs to get out of the collision 
complex, climbing the potential Va{R) until the classical 
turning point Rq where the kinetic energy vanishes, and 
then falling back into the complex, can be approximated 
by twice the time it needs to fall from Rq all the way in. 
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where v{R) = ^2[E -Va{R)]/mr 

RQ^2CQ/mr\/RQ/R^ — 1 is the atom's classical 
velocity as a function of R. The mean collision time At 
follows from averaging over all turning points, starting 
at the equilibrium position R^ of the potential up to 
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The outermost turning point considered, i^max = 
(C6/25^Q™)^/^, is reached if the collision complex scat- 
ters into the energetically lowest closed channel, which 
is the first rotationally excited state of the molecule 
with E = —2Bf^^^^. Explicit values of At for various 
molec ules cover the nanosecond to picosecond regime, cf. 
Using a pure Cq potential as in Eq. (23) instead 



III 



Table 

of our model Lennard- Jones potential and integrating all 
the way to zero instead of stopping at Re is an excellent 
approximation; in the case of Rb + KRb, for example, 
the introduced error is less than 1 %o. 

In cases where P is of order unity or larger (typical 
for heavy molecules), the nuclear spin is almost certain 
to change during the lifetime of the complex. Therefore, 
the nuclear spin degree of freedom also contributes to 
the total DOS. For constructing our statistical models 
we would then use 



p = prv X N^ 



(25) 



where A^nuc denotes the number of nuclear spin states. 
The latter is determined via 



A^n 






(26) 



where /(M, Ml) is the number of possible spin states 
\MAMBfmf) that conserve the total magnetic quantum 
number M for a given Ml. We remark that in prin- 
ciple the ro- vibrational DOS Prv also depends on Ml; 
however, this dependence is negligible such that Eq. (25) 
is a valid approximation. Specific examples of N^uc are 
listed in Table III In calculating A^nuc we include all hy- 



perfine states of the atom but restrict ourselves to sin- 
glet molecular states. Triplet molecular states lead to 
excited electronic quartet and doublet states of the tri- 
atomic collision complex which are connected to the con- 
sidered doublet ground state via avoided crossings [6]. 
The contribution of these states can be assessed by cal- 
culating their DOS at the threshold of the ground state 
doublet potential. Our estimate shows that the overall 
DOS would increase only on the order of 10%. Hence, 
we continue to focus on singlet molecular states solely. 

Our estimate of nuclear spin-changing probability P 
is admittedly crude, and represents at best an order-of- 
magnitude estimate. Nevertheless, all we really need to 
know is whether P is likely to be much larger than unity. 
For heavier molecules, this appears to be the case, and 
we will include nuclear spins in the DOS for our Rb + 
KRb example below. However, for some lighter species, 
such as K + LiK, we expect nuclear spins to be fairly 
well conserved during the collision. 
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FIG. 3. (color online) Elastic partial wave cross sections for 
K + LiK collisions; the top panel shows the corresponding 
cumulative elastic cross sections. The incident channel is the 
absolute ground state; the maximum collision energy is twice 
the van der Waals energy scale ^vdw- 



The determination of the complex lifetime r is not 
influenced by the inclusion of the nuclear spin states 
since both the DOS as well as the number of asymptotic 
channels increase by the same factor: r = 2T:hp/Na = 

27rh{pry X A^nuc)/^^nuc = T = 27rhpry. 



III. HIGHLY RESONANT SCATTERING NEAR 
THRESHOLD 

Using the model described above, we now calculate 
simulated collision cross sections in ultralow energy limit. 



A. Elastic scattering 

We have seen that the density of states can vary widely, 
depending on the particular species we consider. For 
this reason, in this section we will explore two schematic 
cases, where the DOS is either "low" or "high," meaning 
few or many resonant states within the characteristic en- 
ergy scale £^vdw within which the Wigner threshold laws 
hold. 

For the weakly-resonant case, we pick ^^K+^Li^^K. 
For this particular example, we do not expect nuclear 
spin states to be changed (P = 0.2) and hence expect 
only 0.06 s-wave resonances per £^vdw, cf. Tables [ll| and 
|III[ Exemplary elastic partial wave cross sections up to 
L = 4 as a function of the collision energy are shown in 
Fig. [3J the incident channel is the absolute ground state 
of both the atom and the molecule, so only elastic scat- 
tering is possible. As expected from the low DOS, within 
the given energy range (up to 2 x £^vdw) resonances are 



TABLE IV. Density of states p for ^'^Rb +^° K^'^Rb colhsions 
as a function of the partial wave and its magnetic quantum 
number. The total magnetic quantum number is always cho- 
sen such that the absolute ground state is included. 
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encountered only sporadically and are thus well resolved. 
Therefore, the Wigner-law behavior of the elastic cross 
sections is evident: a^^ ex E'^^ for L = 0, 1, and a®^ ex E^ 
for L > 2 [65 . We remark that the resonances found in 
Fig. [3] are determined within our statistical approach and 
hence are representative, not predictive. For a quantita- 
tive description of low-resonant cases such as K + LiK a 
full coupled channel calculation is necessary, at least to 
provide a realistic short-range if-matrix; the long-range 
part then may still be treated by means of MQDT [12 . 
Let us now switch to a high DOS, for which the present 
theory is intended. As a particular example we choose 
^^Rb +^^ K^^Rb collisions [66 for which we expect the 
nuclear spin states to be changed during the form ation o f 
the collision complex (P > 1). Because of Eqs. ( 25|26 ), 
the actual DOS depends on the partial wave considered. 
For higher partial waves L, the 2L + 1 projections Ml of 
the orbital angular momentum allow for a greater variety 
of spin and rotational states that conserve a given total 
magnetic quantum number M. The resulting DOS for Rb 
+ KRb as a function of L, Ml are tabulated in Table [TV| 
Not only does the DOS increase rapidly as a function of 
L, but also all angular momentum projections Ml need 
to be summed to form the final partial wave cross section. 
Since to every Ml a different short-range spectrum is 
attached, this increases the DOS by an additional factor 
of approximately 2L + 1 compared to the case of a single 
Ml. This rapid increase of the number of resonances 
can be observed in Fig. pFb). As indicated by Table IV 



within one £^vdw there are fewer than 10 resonances for 5- 
wave collisions, over one hundred for p-wave, and already 
close to one thousand for d-wave. 

For ultracold temperatures usually all collision pro- 
cesses except for 5-wave are suppressed. In the case of 
highly resonant scattering as we investigate here, how- 
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FIG. 4. (color online) Elastic cross section for Rb + KRb col- 
lisions in the absolute ground state. The maximum collision 
energy is twice the van der Waals energy scale ^vdw- Panel 
(a) shows the sum of the individual partial waves depicted 
in (b). Panel (c) provides a comparison of higher (L > 2) 
partial wave cross sections including [black lines, same as in 
(b)] and omitting (orange/gray lines) the long-range phase 
shift tan^ oc /c^ due to the Cq/R^ van der Waals dispersion 
potential. 



ever, there are plenty of resonance peaks due to higher 
partial waves, cf. Fig. |4ja). Due to the threshold scaling 
of resonance widths [67 , these events are isolated and 
should in principle be well resolvable at the very cold 
end, Ec <C £^vdw- Closer to ^^vdw, on the other hand, 
higher partial waves are not as suppressed and in addi- 
tion resonances start to overlap. Hence, in this regime 
the appearance of the total cross section eventually is 
no longer determined by the background scattering cross 
section with a few resonances on top of it, but rather by 
the interplay of many overlapping resonances. 

The increasing number of resonant states with increas- 
ing partial wave is not the whole story, however. As 
Fig. mb) shows, the number of visible resonances in- 
creases from L = to 1 to 2, but fewer resonances appear 
for L > 2. The reason for this can again be found in the 
Wigner threshold laws. Recall that the elastic scatter- 
ing phase shifts for a van der Waals potential have two 



distinct components. There is a short-range component 
(5sr oc /c^^+^, which vanishes faster with energy for higher 
partial waves because the incident wave function has an 
ever-greater centrifugal barrier to tunnel through. There 
is also a long-range component 5\^ ex k^^ which arises 
from scattering outside the outer classical turning point 
of the centrifugal potential [65 . Since the resonances 
originate in short-range scattering, they appear in (5sr, 
whereby this part of the cross section can be dwarfed by 
(5ir for L > 2. To show this more explicitly, we separate 
out the short-range contribution in Fig. [4Fc) ; within the 
MQDT theory, this amounts to neglecting the e*^ terms 
in Eq. (10). We therefore conclude that, while the num- 



ber of resonances grows rapidly with increasing partial 
waves, nevertheless they are unlikely to be observed in 
the ultracold. 

Instead of the elastic cross sections as a function of col- 
lision energy, presented in Fig. |4J in experimental prac- 
tice one is more likely to measure scattering rate con- 
stants that are thermally averaged. Moreover, often the 
temperature is fixed and, instead, an external magnetic 
field is varied to tune the various scattering channels 
with respect to each other. For these reasons, we pro- 
vide in Fig.l5]the thermally averaged elastic rate constant 
K^i _ ^^ (je^as a function of magnetic field, again for Rb 
+ KRb. Because of the vast difference in energy scales, 
we assume that the short-range physics is independent of 
the applied magnetic field. As a consequence, the colli- 
sion complex probes the short-range resonance spectrum 
with a rate corresponding to the Zeeman shift of the in- 
cident channel, i.e., energies are converted into magnetic 
field strengths via E = — /imag^- For the particular ex- 
ample of ^'^Rb + ^^K^'^Rb with the absolute ground state 
as incident channel, this Zeeman shift is largely deter- 
mined by the magnetic moment of the / = rrif = 1 
ground state of the rubidium atom (/imag = 0-7 MHz/G). 
The resulting densities of states as a function of magnetic 
field are listed in Table |IV| for s-wave collisions as de- 
picted in Fig. [5J it amounts to 5 resonances per Gauss, 
which still should be experimentally resolvable. 

These resonances will naturally wash out with increas- 
ing temperature. In Figure [5] we compare the s-wave elas- 
tic rate constant for three different temperatures, namely, 
1 /iK, 10 /iK, and 100 /iK. As expected, a higher tem- 
perature gradually smoothes the sharp resonance peaks 
found for 1 /iK. Moreover, the few resonances found for 
5- wave scattering as a function of collision energy. Fig. [4j 
have now turned into a dense spectrum of resonances due 
to the large Zeeman shift of the rubidium atom. Since 
sub-/iK temperatures are easily reached in ultracold al- 
kali experiments, we predict that individual resonances 
ought to be observable. 

We remark that in principle Fig. |5] shows also Fano- 
Feshbach resonances that occur within the ro- vibrational 
ground state, i.e., without the need of a highly resonant 
short-range part. However, within the magnetic field 
range shown in Fig. |5| there are only a few such reso- 
nances: the difference between the last and last but one 
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FIG. 5. Elastic s-wave rate constants for Rb + KRb collisions 
in the absolute ground state. Depicted is the thermalized 
rate constant for a temperature of 1 /iK, 10 /iK, and 100 /iK, 
respectively (top to bottom). 



vibrational level of Rb + KRb is on the order of 100 mK. 
The atomic Zeeman shift is 33.5 /iK/G, i.e., within 100 
G it is unlikely to find two different vibrational states of 
the same channel. 



B. Onset of Ericson fluctuations 

Resolving individual scattering resonances requires 
that their mean width f be less than their mean separa- 
tion d = 1/ p. However, as more open channels become 
available, the widths should increase, and when T > d 
the scattering should be in the Ericson regime. For ul- 
tracold collisions, adding more open channels is as sim- 
ple as preparing the molecules in a higher-energy hyper- 
fine state. In this section we therefore explore resonance 
widths as a function of the number Nq of open channels. 

In the limit where the mean resonance width exceeds 
the mean level spacing one might expect the spectrum 
to become increasingly smooth and featureless. Ericson 
showed that surprisingly this is not so [29l|30]. Rather, 
there remains structure that can be probed via the 2- 
point correlation function F{AB) (here written in terms 
of magnetic field strength). 



F{AB) = {a^\B + AB)a^\B)) - {(J^\B)f 



(27) 



where the brackets denote the average over B. In the Er- 
icson regime, this function is predicted to be Lorentzian 



F{AB) (X 



1 + (AB/r)2 • 



(28) 



Since in this section we are only concerned with magnetic 
field dependent cross sections, F has units of magnetic 



field in our case; it easily converts to energy via E = 
— /imag^, though. Similarly, we also consider the DOS to 
be expressed in the magnetic field domain, as in Table 

El 

In nuclear physics, where the number of asymptotic 
channels Na ^ 1, the correlations predicted by Eriscon 
have been nicely demonstrated [ST. Ultracold atom- 
molecule collisions as in the present work, on the other 
hand, are an ideal candidate to investigate the onset of 
Ericson fluctuations for only a few asymptotic channels. 
They possess a large enough DOS for the statistical argu- 
ments to be valid and - most importantly - the number of 
relevant asymptotic channels can be very precisely set by 
choosing the initial hyperfine states of the colliding par- 
ticles. Unlike the treatment in nuclear physics where all 
asymptotic channels are considered open, i.e.. No = Na, 
in ultracold collisions only Nq of them remain open at 
infinite separation of the particles and provide a finite 
outbound flux as i^ ^ oo. Therefore, in the context of ul- 
tracold collisions, we are concerned with No rather than 
Na. Moreover, by setting the short-range coupling pa- 
rameter i?i — 1, we might expect the collision complex 
to decay with the RRKM rate F = Frrkm = No/^Trhp, 
cf. Eq. (fTsl), where p is density of states per magnetic 
field interval. 

In Figure [6] we show the mean widths of resonances as 
a function of the number of open channels. The value of 
the width is extracted from our scattering data as follows. 
We consider magnetic field dependent elastic s-wave cross 
sections for Rb + KRb collisions at a fixed collision en- 
ergy of 100 nK over the range < 5 < lOOG. The 
number of open channels No is set by varying the incom- 
ing channel from the absolute ground state of the system 
(/ = m/ = 1, Mk = -4, MRb = 3/2) to the highest 
possible spin state within the ro- vibrational ground state 
{f = rrif = 2, Mk = -2, MRb = -3/2). We restrict 
ourselves to the conserved total magnetic quantum num- 
ber M = —3/2 that contains the absolute ground state, 
which yields in total 31 possible spin states. After gen- 
erating a model short-range resonance spectrum Ejj^ and 
width matrix W^a^ we compute for each of the No chan- 
nels the magnetic field dependent elastic cross section. 
The mean width of the resonances encountered in these 
cross sections is extracted via their 2-point correlation 
function Eq. (J^TJ) by fitting to the form Eq. ([28|. This 
procedure we repeat for 300 randomly sampled short- 
range resonance spectra and width matrices, for each No. 
To accommodate for the different DOS seen by differ- 
ent hyperfine states, we employ p = |m/|p(m/ = 1). A 
caveat is that, for states with m/ = 0, there is a very 
weak dependence of the incident threshold on magnetic 
field. In this case, only few resonances are seen in the 
magnetic field range < B < 100 G, and we cannot ex- 
tract a width. Thus, in Fig. |6]no data points are shown 
for No = 5 -8 and 21 - 24. 

Figure [6] presents the results of this simulation, show- 
ing the width extracted from Eq. ( 28 ) as a function of the 



number of open channels. The dots represent the mean 
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FIG. 6. (color online) (a) Onset of Ericson fluctuations in 
ultracold Rb + KRb s-wave collisions at a collision energy 
of 100 nK. Shown is the mean resonance width (dots/solid 
line) of elastic cross sections as a function of the number of 
open channels. See text for further details. Insets: Exemplary 
elastic cross sections for No = 2 and No = 9 open channels, 
respectively, (b) Atomic Zeeman shift of the scattering chan- 
nels. The molecular nuclear spin gives rise to an additional 
magnetic field dependence which splits each atomic line into 
a number of sublevels; on the given scale, these sublevels are 
not visible. The channel indices equal the number of open 
channels No. 



value from 300 trials, while the shaded region indicates 
the 1-cr scatter among the trials. The widths are scaled 
by 27rp such that the RRKM width Eq. ^ is equal to the 
number of open channels, i.e., it results in a unit slope 
as a function of Nq (dashed line). The width observed 
in Fig. [6] is unquestionably an increasing function of No^ 
but grows at a different rate for No < 4: and for No >9. 
Focusing on A^o < 4 in Fig. [6) our data shows a slope 
significantly smaller than the RRKM value. This can be 
explained as follows. Assuming that Eq. (fTj) is applica- 



ble, resonances begin to overlap when A^^ > 27r. Hence, 
for No < 4: we are in the regime of isolated resonances 
where T is determined by individual widths and not by 
the collective behavior of many overlapping resonances. 
Nevertheless, we still find a linear behavior, F = ^yNo 
whose slope 7 is determined by the Wigner threshold 
laws. 

The widths of non-overlapping resonances can in prin- 
ciple be extracted from the short-range iC- matrix K^^ . 
More precisely, Eq. (14) yields a mean resonance width 
in the absence of threshold effects of f = 2/7rp, which is 
defined at Rm assuming that particles can freely propa- 
gate beyond this point. Threshold effects, which narrow 
this width, are accounted for within our MQDT treat- 
ment. For a single open channel, and neglecting the po- 
tential resonant infiuence of closed asymptotic channels, 
the MQDT transformation ([9| turns into a simple alge- 
braic equation. Employing further the zero energy limits 
77 ^ and G -^ ( — 1)^+^ [43 , the elastic cross section 
reads 



47r 



a 



A{k,Lf 



k^[l^{-l)L+^K-^]^^A{k,Lf' 



(29) 



For an isolated resonance in the short-range i^-matrix, 
K^r = -{f/2)/{E - E,es)^ Eq. ^ yields a Lorentzian 
shaped resonance in the elastic cross section with a width 
of 

r(7V, = 1) = 7 = A{k, L) X r = ?^iiM) . (30) 

TTp 

The low energy behavior of the MQDT parameter 
A{k,L) is known analytically [431168], 



A{k,L) 



1/2 _ 



R: 



x+1/2 

vdW 



r( 



'.L+1/2 



v^2i-V2(2L + l)!! 



(31) 



^vdw being the van der Waals length i^vdw = 
{2mrCQ/Ji^y^^. For the parameters of Fig. pi (s-wave, 
E = 100 nK), we thus find 27rp x 7 = 0.18. This line 
is shown as a dotted line in Fig. [6] for A^o < 4; it agrees 
quite well with our numerical result. Hence, even in the 
limiting case of isolated resonances. No < 27r, we find 
that the RRKM assumption is true: the decay rate and 
therefore the width of the resonances scales with the in- 
verse of the DOS and is proportional to the number of 
available outgoing channels. 

For No > 9, the resonances nominally overlap since 
T/d = No/2ti > 1 in the RRKM formula. If this in- 
deed places us in the Ericson regime, then the width 
F extracted from Eq. (28) should scale according to 
27rpr = A^o, i.e., should form a line of unit slope in Fig. [6] 
And indeed this is true, apart from an offset, as seen in 
the figure (dotted line for No>9). The source of this off- 
set originates in the Dyson- Wigner distribution of level 
spacings. Qualitatively, the Dison- Wigner distribution 
discourages levels from being close together. The on- 
set of overlapping resonances is therefore deferred until 
higher No. 
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FIG. 7. (color online) Influence of the Dyson- Wigner distribu- 
tion on the correlation function Eq. (27). Shown is the calcu- 
lated width Tout of the correlation function as a function o f th e 
input width Fin for the Lorentzian model spectrum Eq. ( 32 ) . 



Orange dotted line: Resonances are distributed according to 
the Dyson- Wigner distribution Eq. (15). Green solid line: 
The same resonances are normally distributed. The shaded 
areas indicate the standard deviation for the given sample of 
30 individual runs. 



More quantitatively, it can be understood by employ- 
ing a simple model. Instead of the actual cross section, 
we calculate the width Fout of the correlation function of 
a model spectrum that consists of identical Lorentzian 
resonances of width Fin at different magnetic field val- 
ues. 



^mod 



(S) 



Y.\ 



rin/2 



7r(B-B^)2 + (ri„/2)2' 



(32) 



whose locations 5^ are distributed according to Eq. (15). 
We used the same DOS as in Fig. |6] The result is shown 
in Fig. [t] (orange dashed line) as a function of the input 
width Fin of the Lorentzians. For a single resonance, our 
procedure yields Font = Tin- For non-overlapping reso- 
nances. Fin "C 1/p, this can be seen in Fig. [t] as well. 
As the resonances start to overlap, however. Font starts 
to deviate from its linear behavior by showing a smaller 
width than the input. For a strong overlap, the linear be- 
havior is recovered again, but now with a constant offset 
from unity. This is the same qualitative the behavior seen 
in Fig. [6J By contrast, if we repeat the model calcula- 
tion but now for normally distributed resonances instead 
of Dyson- Wigner (solid green line in Fig. [t]), the result 
changes quite drastically: the width calculated from the 
correlation function reproduces the input width. Hence, 
we attribute the offset of the calculated widths to the 



particular statistical properties of the Dyson- Wigner dis- 
tributed resonances. 

We remark that the observation of an autocorrelation 
function with a Lorentzian shape is a necessary but not 
a sufficient condition for the existence of Ericson fluctu- 
ations. As shown, for example, in the case of Helium 
photoionization, one finds a Lorentzian shape of the au- 
tocorrelation function [69 . Looking at the microscopic 
processes in more detail, however, this cross section is ac- 
tually dominated by only a few, well-defined resonances 
not in the Ericson regime. The large number of remaining 
resonances, that would promote the DOS into the Eric- 
son regime, are of too small intensity to contribute to the 
cross section. In our case, on the other hand, there is no 
hierarchy in the intensities and all resonances contribute 
similarly, which makes the autocorrelation function again 
a good indicator for Ericson fluctuations. 



IV. CONNECTION WITH EXPERIMENTS 

Our present model of highly resonant scattering makes 
various assumptions and is by no means meant as a quan- 
titative description. As a qualitative guide, however, it 
does predict various trends from one molecule to another, 
and from one internal hyperfine state to another. In this 
section we summarize what might be gleaned from ex- 
perimental data as they become available. 



A. Density of States 

The first and most obvious measurement would be that 
of the density of states itself. This is most easily mea- 
sured, presumably, by the magnetic field variation of a 
cross section, as in Fig. [5J Even in the highly resonant 
case of Rb + KRb that we explored in detail, we still 
anticipate s-wave resonances spaced an average of ~ 0.1 
Gauss apart, along with perhaps a smattering of higher- 
partial wave resonances at sufficiently low temperature. 
In fact, Eq. (30) tells us that for one open channel p-wave 
resonances at £^ = 100 nK only possess a mean width of 
f = 4 X 10~^ G, which is below typical experimental res- 
olution. For s-wave resonances, on the other hand, the 
required resolution seems experimentally reasonable. For 
molecules in their ground hyperfine state, the most likely 
observable would be loss due to three-body processes near 
each resonance. In this case, one should consider the ef- 
fect on observables due to the width of the three-body 
process, a task we have not attempted here. Otherwise, 
measurements of two-body loss versus field for molecules 
in their first excited state would supply a reasonable ob- 
servable. 

Our estimates of the DOS of various collision part- 
ners rely heavily on the assumption that the entire phase 
space allowed by conservation of energy and angular mo- 
mentum is in fact explored by the collision complex. We 
have argued above that the DOS is surprisingly weakly 
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dependent on the maximum number of vibrational or ro- 
tational states populated (see Table [l|. Nevertheless, our 
estimates of both the DOS and the lifetime of the com- 
plex are almost certainly upper limits. There exists, how- 
ever, already one experiment that constrains the DOS, 
namely, collisions of Rb with ground-hyperfine state KRb 
at sub-/iK temperatures [66 . To infer a DOS from this 
experiment, we assume that there is a universal Rb + 
KRb collision rate given by the quantum threshold model 
of Ref. [TO:, if^b'+KRb = ^(2?i'C'6/m3j,,^j,j,,)i/4. since 
Rb + KRb collisions are stable against reactive losses, the 
overall loss rate is given by the above collision rate time 
the probability that, during the complex lifetime r, an- 
other Rb atom hits the complex and destroys it (the Lin- 
demann mechanism) [22^ . This probability is r times the 

Rb + KRb2 rate ilTgb'+KRb. = A'^^^C./ml^^^^^^^f'^ 
As a result, the loss rate is a quadratic function of the 
Rb density n(Rb), 



function, as in Eqs. (127, 28). Determining this V is pre- 



QT 



-QT 



r(KRb) = rn(Rb)^if;^,VKRb^^b+KRb. 



rn(Rb)^ x 10^^ cm7s 



(33) 
(34) 



The JILA experiment on KRb emphatically did not mea- 
sure a quadratic dependence on the loss rate on n{Rb). 
Nevertheless, we can extract an order of magnitude es- 
timate. From Ref. [66 we infer for a Rb density of 
n{Rb) ^ 0.6 X 10^2 cm-^ a decay rate of r(KRb) ^ 
20 s~^. This sets a rough upper limit to the complex' life- 
time of r ?^ 10 ms. What does this mean for the density 
of states? Suppose the lifetime r is set by the RRKM ex- 
pression Eq. (fTl), assuming only one possible exit channel 
since both the atom and the molecule are in their ground 
states. The upper limit on the density of states is then 
p ^ 10^ //iK. This is a very high density of states, much 
larger than the ~ l//iK we estimate for this case. Hence, 
the experiment at least does not contradict our thinking, 
though it is by no means a measurement of DOS. 



B. Resonance widths 

Within the statistical picture we have outlined, there 
are patterns in the resonance widths as well as in their 
distribution. A main result, which we believe to be quite 
general, is the one in Fig. [6J Namely, the width F de- 
duced from the two-body correlation function of a spec- 
trum grows linearly as a function of the number No of 
open hyperfine channels. The rate of growth is, however, 
different for No < 27r and No > 27r. In the former case, 
dT/dNo is small and is governed by the Wigner threshold 
law. In this case the widths of the individual resonances 
may be difficult to extract; observed widths of very nar- 
row resonances are typically set by the temperature, not 
the intrinsic magnetic-field width of the resonance. 

Vice- versa, in the limit No > 27r, the resonances are 
not individually resolved anyway. Here the spectrum 
becomes a varying background of Ericson fiuct nations, 
characterized by the width F of the two-point correlation 



sumably an easier task, experimentally, than locating and 
measuring the widths of individual narrow resonances. In 
such a case, the simplest version of our model predicts 
that the DOS can be extracted from the slope of F versus 
the number of open channels, via dV/dNo = (27rp)~^. 

In cases where F and the DOS p are both measured in- 
dependently, more information can be determined about 
the microscopic scattering system. Then, as in Fig.[6J the 
dimensionless quantity 27rpF can be directly evaluated as 
a function of No- Our basic model predicts that this re- 
lation will be linear with unit slope, but shifted so as to 
intercept the No axis at some positive value. This shift is, 
as we have argued, a consequence of the Wigner-Dyson 
statistics of level spacings. If, on the other hand, the 
empirical plot of 27rpF versus No intercepts the origin, 
the resonances are more likely normally distributed (see 
Fig. ul. Such data would therefore provide evidence that 
the underlying classical dynamics of the collision complex 
is not chaotic. 

Furthermore, the experimentally determined 27rpF ver- 
sus No may not have unit slope, but may grow more 
slowly, perhaps even nonlinearly. This would indicate 
that the widths are narrower than in the "maximal cou- 
pling" limit we have assumed, in which the coupling pa- 
rameter Ra = I for all asymptotic channels. Determin- 
ing realistic values of coupling constants Ri ^ (or equiva- 
lently, transmission probabilities Ta) would further con- 
strain models of how the ro- vibrational ground state cou- 
ples to the resonant complex, in ways that remain to be 
explored. 



C. Likelihood of Changing Spins 

Estimates of the spin-changing probability P are ar- 
guably the weakest point of our discussion above. Here 
again, experiments should shed light on the true situ- 
ation. For example, if molecules are prepared in their 
second-lowest hyperfine state, one could simply directly 
measure the rate at which molecules are produced in the 
lowest state. When the probability P is of order unity 
(or higher, in our estimates) the spin-changing rate con- 
stant should be on the order of the universal rate con- 
stant i^^^, as described above. However, when P <C 1, 
as we estimate for Na + LiNa or K -h LiK, then the spin- 
changing rate is likely smaller by a factor of P. Even 
an approximate determination of P in this way would 
be a strong and useful constraint on the time scales that 
govern the spin-changing dynamics. 

The argument for whether or not the spin changes de- 
pends on the DOS. Thus we can also infer useful informa- 
tion from a direct measurement of the DOS itself. Within 
our order-of-magnitude estimates we expect the density 
of ro- vibrational states to be nearly universal, i.e., to be 
something like 0.1 5- wave resonances per characteristic 
van der Waals energy ^vdw- For heavier molecules with 
P ^ 1, we expect the density of states to be augmented 
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by the number of spin states, as in Eq. (25). Thus for 



example if the DOS were measured directly for both K 
+ LiK (where spin does not change) and for Rb + KRb 
(where spin can change) , both expressed in units of res- 
onances per £^vdW7 then their ratio would be 



p(Rb + KRb) 
p(K + LiK) 



A^nuc(Rb + KRb) = 31. 



(35) 



On the other hand, if it turns out that the spin cannot 
change in Rb+KRb either, this ratio is closer to unity. If 
instead both collision freely allow spins to change, then 
the ratio would be closer to A^nuc(R'b + KRb)/A^nuc(K + 
LiK) ^ 3.4. Checking this kind of scaling would provide 
information for either verifying or refuting our assump- 
tions about whether and how the spin can change. 



V. CONCLUSIONS AND OUTLOOK 

In the present work we have formulated a theory for 
cold and ultracold atom-molecule collisions that incorpo- 
rates the ro-vibrational Fano-Feshbach resonances in a 
statistical manner while treating the long-range physics 
exactly within MQDT. We provided estimates for the 
densities of states encountered for all non-reactive col- 
lisions involving alkali atoms and heteronuclear alkali 
dimers. The question if during the collision the hyperfine 
states of the collision partners are allowed to change is 
answered by means of a semiclassical approach that es- 
timates the lifetime of the collision complex. As it turns 
out, we expect all systems except the very lightest ones 
to allow for such transitions. This has a great influence 
on the scattering process itself since more resonances are 
accessible, pushing the cross sections at ultracold tem- 
peratures over the limit from showing few to many res- 
onances. Exemplary elastic cross sections as a function 
of the collision energy are provided for K + LiK (no hy- 
perfine change) and for Rb + KRb (change in hyperfine 



sublevels allowed). For the latter, also thermalized rates 
as a function of magnetic field strength for fixed temper- 
ature are shown. Since we assume that the short-range 
physics - and therefore also the ro-vibrational resonances 
- is independent of external fields, the density of states 
becomes a function of magnetic field that is probed by a 
rate corresponding to the atomic magnetic moment. This 
translates the density of states from just a few within the 
Wigner threshold limit to many per Gauss. 

One of the intriguing aspects of the considered ul- 
tracold collisions is that the initial states can be very 
well controlled. This allows one to tune the number of 
open channels very precisely and opens the opportunity 
to probe the onset of Ericson fiuctuations. The latter 
are well-known in nuclear physics where one encounters 
usually a large number of (open) asymptotic channels. 
Here we showed that the scaling law predicted by Ericson 
should be nicely observable in ultracold atom-molecule 
collisions. Moreover, by limiting the number of open 
channels to a small number, one can switch between be- 
ing in the Ericson regime and the regime of isolated res- 
onances. 

The Ericson fiuctuations only depend on the density of 
states and the number of open channels. Since the latter 
are fixed by the choice of the initial state, a measurement 
of the Ericson fiuctuations should in principle allow for 
an experimental determination of the density of states - 
an intriguing possibility which is subject of future inves- 
tigations. Also, in the present work only elastic processes 
are investigated. The extension to inelastic or chemically 
reactive ones is straightforward and promises further in- 
sights in the physics of highly resonant scattering. 
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